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ABSTRACT 

Guidance for the formulation of robust, multiaxial, constitutive models for advanced materials 
is provided by addressing theoretical and experimental issues using micromechanics. The 
multiaxial response of metal matrix composites, depicted in terms of macro flow/damage 
surfaces, is predicted at room and elevated temperatures using an analytical micromechanical 
model that includes viscoplastic matrix response as well as fiber-matrix debonding. Macro 
flow/damage surfaces (i.e., debonding envelopes, matrix threshold surfaces, macro “yield” 
surfaces, surfaces of constant inelastic strain rate, and surfaces of constant dissipation rate) are 
determined for silicon carbide/titanium in three stress spaces. Residual stresses are shown to 
offset the centers of the flow/damage surfaces from the origin and their shape is significantly 
altered by debonding. The results indicate which type of flow/damage surfaces should be 
characterized and what loadings applied to provide the most meaningful experimental data for 
guiding theoretical model development and verification. 

INTRODUCTION/BACKGROUND 

Aerospace structural components are commonly subjected to thermomechanical service 
conditions, involving multiple temperature levels, thermal transients and mechanical loads, 
severe enough to result in significant nonlinear (e.g., inelastic or damage induced) deformation. 
The accurate prediction of the stress, strain and temperature fields within these structural 
components depends strongly upon accurate mathematical representations (constitutive 
equations) of the nonlinear behavior of structural materials (e.g., monolithic alloys or advanced 
composite systems) at high temperature. To be generally applicable, these constitutive equations 
must be expressed in multiaxial form and be appropriate for all anticipated modes of mechanical 
and thermal loading (e.g., cyclic, non isothermal, nonproportional, etc.). This is an imposing task, 
nevertheless significant advancement has been achieved over the past two decades with the 



advent of unified viscoplastic constitutive models (e.g.. Walker (1981), Bodner (1987), Helling 
and Miller (1987), Robinson and Duffy (1990), Freed and Walker (1995), Arnold et al. 
(1996a,b)). Historically, with the exception of Spencer (1972), uniaxial constitutive response has 
been the focal point of research on anisotropic materials. Recently, attention has being focused 
on describing the multiaxial constitutive behavior (e.g., Robinson and Duffy (1990), and Saleeb 
and Wilt (1993)), for new generations of advanced structural materials that possess strong 
directional characteristics (e.g., directionally solidified metals, metal and intermetallic matrix 
composites as well as ceramic and polymeric composites). Such advanced materials offer the 
designer the additional advantage of tailoring his design to have increased strength, stiffness, 
conductivity and the like in one or more preferred directions. Under high temperature 
thermomechanical service conditions these advanced materials still exhibit all the time- 
dependent hereditary behavior of more conventional alloys but now their directional nature adds 
yet another level of complexity. Consequently, the further development and validation of 
multiaxial viscoplastic models for metal matrix composites is essential if these advanced 
materials are to be utilized by the aerospace and power generation industries. 

A key ingredient in the formulation of robust, multiaxial, constitutive models is the close 
interaction between experimentalists and theoreticians in establishing the required developmental 
program. Three types of experimentation (Robinson (1985)) are necessary to support the rational 
formulation of constitutive models for high temperature structural material systems (be they 
monolithic or composite). These are 1 ) Exploratory tests that guide the development and examine 
the fundamental concepts and mathematical constructs embedded in the theoretical framework, 
2) Characterization tests that provide the required data base for determining the specific 
functional forms and material parameters to represent a particular material over a given range of 
conditions, and 3) Verification tests , (often prototypical in nature) that provide the ultimate test 
of a constitutive model through comparison of actual structural component response with 
analytical predictions based on a given model. As results from these tests should ideally provide 
feedback for subsequent developmental efforts, the best verification tests may be those of simple 
structures such as beams, plates, shells and bar structures under prototypical conditions of 
temperature, stress, strain-rate, etc. 

The present work will primarily focus upon practical and theoretical issues relative to 
conducting and analyzing tests of the exploration type, i.e., determination of multiaxial 
flow/damage surfaces for titanium matrix composites (TMCs). However, it is important to 
realize that such exploratory tests on composite materials (which are actually structures) can in 
themselves be viewed as verification tests (for both the constitutive model and homogenization 
(micromechanics) technique), provided the constitutive equations being verified are developed 
utilizing only constituent material properties. This distinction will be developed further within 
the body of the report. The necessity of conducting such fundamental multiaxial flow/damage 
surface tests can be best understood in the context of classical plasticity; in as much as the 
concept of a yield surface plays a central role in that it delineates the elastic and inelastic 
deformation regions of the material. The existence and description of this yield surface is 
precisely that which allows a consistent multiaxial statement of plastic flow (flow law) to be 
written. Mathematically, this is done by making use of the fundamental assumption that the yield 
surface has the properties of a potential (normality), e.g. 
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and thus leads to an associated form of the flow law through differentiation of /, that is 
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where, o t] , a ip £ l ip T, and k are the applied Cauchy stress, back (internal) stress and inelastic 

strain tensors, temperature and yield stress, respectively. 

Experimentally, yield surfaces are typically mapped out using an equivalent inelastic strain 
definition and a rather arbitrary target value of permanent set. The impact of the magnitude of 
this target value will be indicated subsequently. Multiaxial testing alone provides a quantitative 
description of/, thereby allowing a correct mathematical form of the flow law to be deduced. To 
date significant testing has been done in support of the concept of a yield surface and the 
normality condition, albeit at relatively low temperatures and for monolithic materials (e.g. 
Phillips (1974), Liu and Greenstreet (1976), Ellis et al. (1978) and Oytana et al. (1982)). Note, as 
we are concerned here with primarily examining composite materials, “yield” will be defined to 
have occurred whenever deviation from proportionality (DFP) has been exceeded by a specified 
target value. In this way the inelastic response can encompass both matrix inelasticity and 
internal damage resulting from fiber-matrix debonding. 

At elevated temperatures, metallic alloys and TMCs typically exhibit strong time-dependent 
deformation behavior and the concept of a yield surface, in the classical sense, generally breaks 
down as we move into the realm of viscoplasticity. However, analogous geometrical and 
thermodynamically based concepts such as surfaces of constant inelastic strain rate (SCISRs) and 
surfaces of constant dissipation rate (SCDRs) 1 can and have been postulated to play the same 
central role in viscoplastic constitutive theories as yield surfaces do in classical plasticity. In fact, 
given a complete (i.e., fully associative) potential based framework, as that proposed by Arnold 
and Saleeb (1994), the experimental multiaxial identification of the Gibbs and dissipation 
potential would completely define the ensuing compliance operator and flow and evolution laws 
for the inelastic strain and internal state variables, respectively. This fundamental approach to 
the formulation of a consistent multiaxial theory is in stark contrast to the more typical ad hoc 
and generally inadequate approach of extending uniaxial constitutive theories for multiaxial 
conditions by simply placing bars over the pertinent variables and declaring them “effective” 
values. This systematic approach becomes even more important when one deals with materials 
that have preferred directions (e.g. TMCs), as the anisotropy and potential for internal damage 
(i.e., since debonding is prevalent in TMCs, where it can occur prior to plastic flow in the matrix) 
further complicates the shapes of these surfaces in various stress spaces. Many researchers have 
documented the effects of debonding on the overall deformation response of TMCs (e.g., 
Johnson et al. (1990), Jansson et al. (1991), Majumdar and Newaz (1992), Lerch and Saltsman 
(1993), Lissenden et al. (1996a)). 

Both of the above definitions (i.e., for SICSRs see Battiste and Ball (1986) and for SCDRs 
see Clinard and Lacombe (1988)) have been used to describe the form of the dissipation potential 
at elevated temperatures for monolithic materials. However, the utility (both theoretically and 
experimentally) of using a SCISR or SCDR concept to define the requisite state dependence of 


1 SCDRs are defined to be 2,// where I y is the effective deviatoric stress tensor (i.e., 2, ; = s,j - a :J , with s u being the 
applied (external) deviatoric stress and the back (internal) deviatoric stress) and ej y the inelastic strain rate tensor; and 
SCISRs are defined as £ y £ y . 
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the dissipation potential requires careful consideration as these “definitions” may or may not be 
equivalent, particularly with regard to quantifying the multiaxial behavior of anisotropic 
materials. To substantiate this statement, let us first consider an isotropic, monolithic, material 
having a dissipation potential 0(F,G) that depends on the external stress through the scalar 
function F and on the internal variables (the back stress in this case) through the scalar function 
G (Robinson and Ellis (1986) and Hopkins (1990)). From normality, the flow law becomes, 

. i _ dO _ dO dF 
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where J 2 ==-^-ZyX,y and 7 3 =^J. i j'Lj k L ki are the second and third invariants of the effective 

deviatoric stress, respectively, c is a constant, and f (F) is an unspecified function, then for a 
fixed inelastic state , i.e. G=constant » 0, 
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where J A = ^'L ik 'L k j'L il 'L[j . Thus SCDRs and SCISRs are both surfaces of constant Q(F,G), i.e. 

F = constant, if and only if the material is a purely J 2 material, i.e. c = 0, and then the SCISR 
becomes 
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Consequently, because SCDRs are independent of the mathematical constructs of F, they are 
the most theoretically meaningful measurements when desiring to quantify the state dependency 
of the dissipation and assess the applicability of the important normality condition (i.e., 
convexity of this potential). 

Next, consider a macromechanical viscoplastic model that describes the inelastic behavior of a 
unidirectional metal matrix composite which possesses transverse isotropy (Robinson and Duffy, 
1990). Here the scalar function F is defined as 
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where the following definitions apply: 
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Dij is formed from the self-product of the unit vector pointing in the fiber direction, and r| and co 
denote the threshold stress ratios for shear and normal loading, respectively. The subscript L 
refers to the longitudinal direction and T the transverse direction. Similar to the isotropic 
material case, it can be shown that SCDRs are proportional to surfaces of constant F (i.e., 
£2=constant) whereas SCISRs are not. In fact, even if we restrict the problem to a longitudinally 
reinforced tube subjected to only axial-torsional loading, the SCISRs are not necessarily 
proportional to surfaces of constant F. SCDRs and SCISRs for this case are, 


Z ij £‘ lJ =2f\F)4(F + \) 
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respectively, where cr denotes the axial stress component and x the shear stress component, 
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and F reduces to 
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The fact that SCISRs are not necessarily proportional to surfaces having a constant value of F 
is significant because the desirable convexity and normality features associated with constant F 
(i.e., £2=constant) surfaces are not guaranteed for these SCISRs. Experimentally however, 
SCISRs have been thought to possess a distinct advantage over SCDRs for TMCs, as they can be 
determined from multiple probing experiments on a single specimen because only a negligibly 
small change in material state occurs during each probe. As will be discussed in more detail 
later, this may not be the case for SCDRs, where due to the anisotropy, a large inelastic strain rate 
must accompany the relatively low stress in the weak direction. Furthermore, stress 
concentrations leading to very localized inelastic deformation could conceivably preclude multi- 
probe SCDR or SCISR experiments on a single specimen. 

With this background in mind, our objective in this paper is to analytically (using the method 
of cells micromechanics model) generate macro flow/damage surfaces for a representative 
titanium matrix composite system (i.e., SCS-6/TIMETAL 21S) and critique the resulting 
geometric shapes (given various definitions of these surfaces). This critique will be conducted 
relative to the theoretical implications when desiring to represent such macro surfaces using a 
macromechanics based constitutive model and the experimental implications when attempting to 
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experimentally measure them. It is hoped and anticipated that this work will provide insight and 
guidance for future experimental programs aimed at developing macro anisotropic constitutive 
models as well as multiaxially validating currently proposed micromechanical theories. 

MICROMECHANICAL MODEL 

The method of cells micromechanical model (detailed in the monograph by Aboudi, 1991) 
was used by Pindera and Aboudi (1988) to predict yield surfaces for MMCs. Micromechanics 
was also employed by Dvorak (1991) to determine yield surfaces for MMCs using a bimodal 
plasticity theory (Dvorak, 1987). In both of these efforts perfect bonding was assumed between 
the fiber and matrix. Aboudi (1989) incorporated debonding due to shear loading into the 
method of cells and predicted initial yield surfaces and debonding envelopes. The current work 
employs a more general debonding model and employs different concepts to define flow/damage 
surfaces. 

The constitutive relations and/or parameters used to describe the fiber, matrix, fiber-matrix 
interface, and the overall composite are presented here. A model TMC system (SCS- 
6/TIMETAL 2 IS), having a unidirectional continuous silicon carbide (SiC) fiber reinforcement, 
is used throughout this paper. Additionally, given that a model material concept (Arnold and 
Castelli (1995)) holds, the same approach and conclusions found here should apply to other 
TMCs systems as well. 

Silicon Carbide Fibers 

Textron’s high-strength, high-stiffness, continuous SiC fibers (i.e., SCS-6) are approximated 
as being isotropic and linear elastic. Furthermore, since the variation of elastic properties with 
temperature is relatively small, the following isothermal elastic properties are used throughout: 
Young’s modulus, E = 400 GPa; Poisson ratio, v = 0.25; and coefficient of thermal expansion, 
CTE = 4.5x1 O' 6 °C 1 . 

Titanium Matrix 

The titanium alloy TIMET AL 2 IS® (Ti-21s) is a metastable beta strip alloy, containing 
approximately 21% alloying additions, that has high strength as well as good creep and oxidation 
resistance. Ti-2 Is is commonly used in advanced metal matrix composites, hence its (isotropic) 
viscoplastic response has been characterized for the model of Bodner and coworkers (Chan et al., 
1988), Chan and Linholm, 1990)) by Kroupa and Neu (1993) as well as a generalized 
viscoplasticity with potential structure (GVIPS) model (Arnold and Saleeb, 1994) by Arnold et 
al. (1996a,b). The GVIPS model of Arnold et al. (1996a, b) has been selected for the current 
study. The model has one tensorial internal variable, the back stress, whose evolution includes 
nonlinear kinematic hardening and both thermal and strain induced recovery mechanisms. The 
use of a nonlinear compliance operator in the evolution law permits the flow and evolution laws 
to be fully associated, which greatly influences the multiaxial response under nonproportional 
loading (Arnold et al., (1995)). 

The flow and evolution equations are, respectively, 


6 



if F<0 


•U 3 s!, I, 


f=r » J 

2p 2 

Qijmn^mnkl^kl 


if F >0 


if flyly < 0 
ifflyZy > 0 


&ij - L ijki{\i ~ ®kiTy 

where the dot represents time differentiation, A,j is the deviatoric back strain tensor, a i} the 
deviatoric back stress tensor, T is temperature, and 
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G = -4- (back stress function) 
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I 2 =—a i ja i j, J 2 = — XyZy (stress invariants) 

where £,,*/ are elastic stiffness coefficients, H v [.] is the Heaviside unit step function, and <•> are 
Macauley brackets. The temperature-independent material parameters are: Ko, n, B\, p and q, 
while the temperature-dependent material parameters are: K, p o , Bo, R a , and (3. Interpolation 
functions defined by Arnold, et al. (1996a)) are employed to determine the material parameters 
for Ti-21s at temperatures other than the reference temperature, 650°C. A limitation of the 
material characterization is that above 704°C, material parameters are taken to be those at 704°C. 
Material parameters at various temperatures are presented in Table 1 . 

Fiber-Matrix Interface 

As alluded to in the introduction, fiber-matrix debonding significantly affects the stiffness and 
strength of TMCs. The weak interfacial bond is due to the multi-layer coating applied to the 
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fiber prior to consolidation. In this section, interfacial constitutive equations between tractions, 
Ti, and displacements, w„ for an interface subjected to combined normal and tangential loading 
are summarized from Lissenden (1996). Experimental results suggest that interfacial debonding 
initiates upon the traction attaining a critical value - defined as the interfacial strength, T 0 . For 
the simple case of normal separation, the interfacial strength and ductility are sufficient to 
characterize interfacial separation. The ductility, u 0 , represents the interfacial displacement at 
which the traction vanishes. If, after initiation, debonding proceeds in a smooth fashion, the 
traction-displacement relation can be taken to be a cubic polynomial (Fig. 1): 



+ 2 




\ U oJ 


( 11 ) 


To generalize eq. (11) for combined normal and tangential loading, an x,n,t interface 
coordinate system is defined having the x-direction parallel to the fiber, n-direction normal to the 
interface, and f-direction tangent to the interface and orthogonal to the fiber direction. The 
traction and displacement components are normalized as follows: 



where a defines the tangential-to-normal ductility ratio and A defines the shear-to-normal 
strength ratio. In addition, a normalized Coulomb frictional traction is defined 

' 03) 

where (i is the coefficient of static friction. A quadratic interaction of the tractions is employed 
to predict the initiation of debonding to occur when 

r 2 2 1 1/2 

r={(i , ,-sgn(r x )J» +{5) +{i;-sgn(^)f / ) } >i 04> 


where 


sgn(x) = 



for x < 0 
for x > 0 


(15) 


and Macauley brackets are used in the normal traction term because compressive loading is not 
considered to be detrimental to the condition of the interface. Macauley brackets are used in the 
tangential traction terms to ensure that tangential interfacial displacements do not occur until 
after the static frictional force has been overcome. 

Once debonding has initiated, the tractions can be related to the displacements by generalizing 
eq. (1 1) and including frictional effects as follows: 


T n = F*(u max )u n , for tensile separation 
T x = F*(u max )u x +sgn(u x )T f 

^ = F *(«ma X )“f + Sgn (u,)T f 


>for 0 < u < 1 


(16) 
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and 


where 


T n = 0, for tensile separation 


T x =sgn(ii Jt )7y 
T, = sgn(i7,)7} 


(for u > 1 
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represents the normalized resultant interfacial displacement and 

F*(u) = (l -3« 2 +2« 3 )y/i7 (19) 

represents the generalized form of eq. (11). The maximum value of u = u max is utilized in eq. 
(19) to prevent the interface from being self-repairing upon unloading. The primary difference 
between this model and that of Needleman (1987) and Tvergaard (1990) is the inclusion of a 
debond initiation criterion. 


Homogenization 

The method of cells (Aboudi, 1991) is a versatile micromechanical model for predicting the 
overall response of composites. Overall elastic properties are determined in closed-form. The 
microstructure is assumed to be doubly periodic, therefore it is sufficient to consider a unit cell 
containing a single fiber. The unit cell consists of one fiber subcell and three matrix subcells, 
although a recent generalization of the method permits an arbitrary number of subcells, enabling 
the analysis of various microstructures (Paley and Aboudi, 1992 and Aboudi, 1995). For the 
current work the standard version of the method of cells is adequate, the generalized method will 
be used in future work to study the effects of microstructure on flow/damage surfaces. Subcells 
are represented by rectangular geometries strictly as a mathematical convenience. The square 
unit cell results in a material having cubic symmetry, no attempt has been made in the current 
work to force the material to be transversely isotropic. Since equilibrium of tractions and 
compatibility of displacements are imposed at the interfaces in an average sense, the actual 
geometry of each subcell is only felt in a weighted average sense, i.e., the “square” fiber is really 
only a pseudo-square fiber as only the ratio of l/h is important. 

A linear displacement field in each subcell is sufficient for determining overall response. The 
displacement field is continuous except at the fiber-matrix interface, where average interfacial 
displacement components are defined. Additionally, tractions must be continuous across subcell 
interfaces. It is recognized that the subcell inelastic strains are approximate, but also that the 
numerical efficiency of the model far exceeds that of other, perhaps more exact methods such as 
finite element analysis using a large number of elements. 

Altogether, homogenization provides 24 equations, but there are 30 unknown variables. The 
fiber-matrix debonding relations presented in the previous section provide the additional six 
equations. McGee and Herakovich (1992) showed how to reduce the full system of equations to 
six with the interfacial displacements as the only unknown variables by considering different 
types of applied loading. 


9 



RESULTS 

Aboudi’s micromechanical model was applied to an SiC/Ti composite having a fiber volume 
fraction of 0.35 and the constituent properties discussed in the previous section. Residual 
stresses associated with processing were accounted for by cooling the composite down from a 
stress-free temperature of 815°C to room temperature in two hours. Subsequently, strain 
controlled probing tests were conducted at a relatively slow rate, lOp/sec, to simulate the rate that 
might be used in an experiment. The material was in a pristine state (as processed) at the 
beginning of each probe and 1 80 probes were employed to determine each flow/damage surface. 
The interface properties are shown in Table 2 and were chosen to be representative of an SiC/Ti 
interface (Lissenden, 1996). However, in some cases the interfacial strength, T 0 , has been 
allowed to take on different values, e.g., 100 MPa and 300 MPa. 

To illustrate the marked difference in the effects of matrix inelasticity and fiber-matrix 
debonding. Fig. 2 shows the transverse tensile response (90°) for; (i) matrix inelasticity/perfect 
bonding (a continuous displacement field), (ii) an elastic matrix/debonding fiber-matrix interface, 
and (iii) matrix inelasticity/ debonding. Residual stresses associated with processing have been 
considered in all three predictions. Clearly, debonding causes a deviation from proportionality 
(DFP) prior to the onset of matrix inelasticity. Note that the unloading modulus (stiffness) has 
been degraded for the two cases in which debonding occurred. The precise representation of this 
unloading response, however, is relatively unimportant in the present work as only the loading 
portion of each probing experiment is utilized. 

The stress-strain responses for axial, transverse, and axial shear loadings are shown in Fig. 3 
for three temperatures of interest: 23, 300, and 500°C. Residual stresses have been considered 
and the (weak) interface properties shown in Table 2 were employed and assumed to be 
temperature-independent. While the tensile axial response is independent of fiber-matrix 
debonding, both the tensile transverse and shear responses are greatly affected by debonding as 
one would expect. In this interface model, as in others, an instability associated with the 
redistribution of stress from the fiber to the matrix as the interface debonds can occur, as it does 
for transverse tension at 300 and 500°C. No significant instability is evident at 23°C, where the 
interface properties were chosen to match experimental results. This suggests that more accurate 
results could be obtained by determining the interface properties as a function of temperature,. 

Flow Surfaces for Perfect Bonding 

Composite, or macro, flow surfaces of constant deviation from proportionality (SCDFPs) in 
transverse-axial ( 022 - 011 ), transverse-transverse (O22-O33), and shear-axial (012-On) stress spaces 
are shown in Fig. 4 for the perfect bond condition and target deviations of 10|i and 0.2% (2000 
|x). The deviation from proportionality is due entirely to matrix inelasticity for perfect fiber- 
matrix bonding. Thus, for this case, SCDFPs are akin to yield surfaces, but in general could also 
be associated with internal damage. Overall SCDFPs are defined by the equivalent inelastic 

strain invariant, 

lOjx definition gives an indication of first yielding and has been shown to be significant for 
experimental yield surface determination (e.g., Phillips et al. (1972), Liu and Greenstreet (1976), 
Helling et al. (1986), Lissenden et al. (1996b)) in monolithics. Likewise, the value of 0.2% is the 
commonly used strain offset definition for monolithic materials and provides an indication of 


= A |— e\j£\j , attaining a critical value; i.e., 10|i and 0.2% as in Fig. 4. The 
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significant plastic deformation. In reality, the 0.2% definition is frequently not applicable to 
reinforced materials because the reinforcing can fail prior to achieving an overall permanent 
strain of 0.2%, particularly for loading in the direction of the reinforcing. This however, 
underscores the importance of selecting an appropriate flow definition. The tensile fiber fracture 
envelope (for a fiber strength of 4000 MPa) is also shown in Fig. 4 to illustrate that 0.2% is not 
realistic for SiC/Ti-21s, nonetheless it is of academic value for comparing localized inelasticity 
with gross overall inelasticity. 

Applying the target values of 10(t and 0.2% to our model TMC system, one can see not only a 
difference in size of the resulting flow (yield) surfaces depicted in Fig. 4, but in the G22-G11 and 
G22-G33 stress spaces a noticeable change in shape is also evident. Table 3 presents a comparison 
between the local (subcell) equivalent inelastic strains and their macroscopic (overall) 
counterparts at various locations on the SCDFPs of Fig. 4. For the selected points, the ratio of 
the maximum local equivalent inelastic strain to the value of the macroscopic SCDFP ranges 
from 4.1 (point E) to 6.1 (point A) and 2.4 (point D) to 3.6 (point H) for 10(1 and 0.2% SCDFPs, 
respectively. The fact that the maximum-local-to-macroscopic equivalent inelastic strain ratio is 
larger for the 10|i target value than for the 0.2% case provides the proper indication that matrix 
inelasticity is more localized for the 10(1 SCDFP. The comers appearing on the 10|i SCDFPs 
(Fig. 4) are associated with the change in location of the localized inelastic strain, i.e., on one 
side of a comer matrix inelasticity occurs in one subcell while on the other side of the comer it 
occurs in a different subcell. These comers are an artifact of the approximate microstructure 
employed by the method of cells and were also observed by Pindera and Aboudi (1988). Future 
work will investigate more accurate microstructural representations. Finite element analysis 
(Dvorak et al., 1973) has validated the overall SCDFPs obtained by the method of cells. 

The implications of these results for an experimental program reside in the fact that if it is 
desired to limit local (constituent) equivalent nonlinear strains (so as to minimize any change in 
internal state) to 50(1, for example, it would be necessary to choose an overall DFP value less 
than 10(1 in the G22-G11 stress space while an overall DFP value slightly greater than 10(l could be 
employed in the G12-G11 stress space. In either case, the in-situ, unmeasurable (using 
extensometry or strain gages) local, nonlinear strains would be approximately 5 times greater 
than their macroscopic (measurable) counterparts. 

Utilizing the alternative SCDR concept to define overall flow surfaces, the predicted SCDRs 
from the current micro-model and surfaces of inconstant (eq. (4)) from the macro-model of 
Robinson and Duffy (1990) are shown in Fig. 5 for the case of no residual stresses. The strengths 
of anisotropy required for the macro-model were obtained by using the micro-model to simulate 
normal and shear loading experiments, yielding a> = 1.7, 17 = 0.83, and k T = 610 MPa. The value 
of F was chosen such that the micro-predicted and macro-predicted overall surfaces agreed in the 
G12-G11 stress space. The two overall SCDRs agree reasonably well in the G22-G11 stress space, 
but less so in the G22-G33 stress space where transverse stresses dominate. This observation is 
supported by Arnold et al. (1993) in related work dealing with creep-plasticity interaction 
wherein it was shown that non-associated models (such as Robinson and Duffy, 1990) can 
accurately represent longitudinal behavior and yet are inadequate for transverse stresses, while 
associated models of the GVIPS class (see Arnold and Saleeb, 1994, and Arnold et al., 1996a,b) 
provide more consistent results. It is not obvious why applying a small amount of G22 causes the 
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dissipation rate target value to be achieved at a larger an (Fig. 5), but this is indeed the case as 
was also observed by Pindera and Aboudi (1988) and Dvorak et al. (1973) for yield surfaces. 

The centers of all surfaces in Fig. 5 are located at the origin because residual stresses have not 
yet been included. If we now account for residual stresses in the micro-model the centers of the 
surfaces shift as shown in Fig. 6. The residual stresses are observed to strengthen the composite 
under equi-biaxial transverse tension, but weaken it under equi-biaxial transverse compression 
(Fig. 6c). In the macro-model we attempt to account for residual stresses by equating the initial 
back stress to the center of the predicted micro-model surface, here an = -241 MPa. The 
comparison shown in Fig. 6 indicates that this is a viable macro-modelling technique to account 
for a material having different flow stresses in tension and compression. 

Applied loading (strain) rate can greatly affect flow behavior, particularly at elevated 
temperatures. The SCDFPs shown in Fig. 7 for loading rates of 2, 20, 200, and 800 p/sec at 
500°C illustrate the rate dependency that one would expect for such a surface (here one should 
not think of a yield surface, but rather a surface of constant inelastic strain). Since inelastic 
deformation takes time to occur, more inelasticity will result from slower loading rates. On the 
other hand, SCDRs (and SCISRS) depend on the inelastic strain rate and thus are affected quite 
differently by loading rate. SCDRs are observed to decrease in size, and change shape, as the 
loading rate is increased. Furthermore, the loading rate is an upper bound on the inelastic strain 
rate and if the loading rate is too slow (e.g., 25 p/sec in Fig. 7b) it may not be possible to 
determine an entire SCDR. As will be discussed subsequently, a discontinuity in the inelastic 
strain rate occurs at the threshold stress due to the formulation of the viscoplastic material model. 
Once the loading rate gets fast enough, this discontinuity will cause the target dissipation rate 
value to be exceeded at the threshold stress, which explains why SCDRs for 1000 |i/sec and 
10000 p/sec are essentially the same and suggests that for extremely fast loading rates SCDRs 
are rate-independent. 

Flow/Damaqe Surfaces 

Damage envelopes for metal matrix composites have previously been predicted by Aboudi 
(1989). Here, the concept of a damage envelope has been combined with that of a yield or flow 
surface resulting in SCDFPs (rate-independent) or SCDRs and SCISRs (rate-dependent). In the 
current work we are defining damage to be fiber-matrix debonding, other forms of damage such 
as matrix cracking, fiber breakage, and void growth are specifically excluded. The effects of 
local events (i.e., debond initiation and subcell stresses attaining their threshold value, F 0) on 
overall SCDFPs are shown in Fig. 8 by superimposing damage initiation (eq. (14)) envelopes and 
matrix threshold (eq. (10)) loci on overall lOp SCDFPs in transverse-axial (a 22 -<Ju), transverse- 
transverse (a 22 -CJ 33 ), and shear-axial (ai 2 -Gn) stress spaces at 23°C. Be aware that the initiation 
of debonding gives no indication of how the interfacial displacements accrue. It is clear from 
Fig. 8, however, that axial loading and compressive biaxial transverse loading do not initiate 
debonding. (As will be shown later, axial compression can lead to debonding associated with 
Poisson expansion.) The lower extreme of the debonding envelope in the <3ii-<5\\ stress space 
(Fig. 8b) is associated with overall (macro) Poisson expansion due to compressive transverse 
applied stresses. Here, for an interfacial strength of 200 MPa, the DFP target value is attained 
prior to the matrix threshold when debonding occurs due to transverse loading but not due to 
shear loading, thereby once again suggesting that debonding affects the transverse tensile 
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response more than it does the shear response, as one might expect. Thus, the overall 
flow/damage SCDFP is essentially bimodal for normal loadings, in that it is the intersection of 
the debonding envelope and matrix threshold. However, in the presence of shear loading, an 
interaction between debonding and matrix inelasticity lead to deviation from proportionality. 

The effect of temperature on 10p SCDFPs is illustrated in Fig. 9 where flow/damage surfaces 
have been predicted for 23, 300, and 500°C. As expected, in all three stress spaces the size of the 
SCDFP decreases with increasing isothermal conditions and, as the interface properties were 
taken to be temperature-independent, the only effect temperature has on the initiation of 
debonding is in the magnitude of residual stresses produced (as these stresses .which serve to 
clamp the interface together, relax as the temperature increases). As depicted in Fig. 3, the 
predicted deviation from proportionality associated with debonding occurs within a relatively 
small range of transverse and shear stresses. 

Similarly, Fig. 10 shows how interfacial strength affects SCDFPs at 23°C. In addition to 
SCDFPs for interfacial strengths of 100, 200, and 300 MPa, the perfect bond condition is also 
shown. In 022-011 stress space (Fig. 10b), debonding due to Poisson expansion affects the 
SCDFP for compressive axial and transverse stresses in the case of an interfacial strength of 100 
MPa, but not for higher interfacial strengths. As the interfacial strength of the composite is 
increased, the size of the SCDFP increases in the region where debonding has an influence with 
the perfect bond condition as the extreme (Fig. 10). 

The effect of interface ductility, u 0 , was also studied, but found to have a relatively small 
effect on the size and shape of the flow/damage surfaces. This was expected, since interface 
ductility only affects the rate at which stresses are redistributed from the fiber to the matrix upon 
debonding. In the <122-0 1 1 and O33-O22 stress spaces, inelasticity is observed almost immediately 
after debond initiation and only an unrealistically high ductility would impact the inelastic 
surfaces. Greater interface ductility did result in slightly larger inelastic surfaces in the On-On 
stress space (not shown). 

Now consider the construction of rate-dependent SCDRs and SCISRs at 500°C (Fig. 1 1 and 
13). Both types of surfaces depend on the overall inelastic (including debonding effects) strain 
rate. Here it was found that the applied strain rate of lOp/sec was too slow, in that it placed an 
upper bound on the inelastic strain rate and severely limited the range of values that could be 
used for complete SCDR and SCISR determinations. Thus, the applied probing strain rate was 
increased to lOOp/sec for all subsequent results. With respect to SCDRs, three target values, i.e., 
1, 20, and 35 Pa/sec, are shown in three stress spaces at 500°C in Fig. 11. As with SCDFPs, 
SCDRs exhibit straight segments where fiber-matrix debonding dominates and fairly sharp 
comers at the intersection of regions dominated by matrix inelasticity and those dominated by 
debonding. 

It is important to note that these fairly geometrically complex macro SCDRs, were predicted 
from a micro-model having relatively simple constituent constitutive relations (i.e., a J 2 
viscoplasticity theory for the matrix and a simple multiaxial interface model). Conversely, the 
required mathematical form of a macro-model would have to be quite sophisticated to obtain the 
predicted SCDR representations shown in Fig. 1 1 . The primary difficulty in modelling TMCs at 
the macroscale would appear to be in properly accounting for strain and stress localizations and 
the differences in tension and compression associated with damage. 
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Even if viable mathematical constructs for the dissipation potential can be shown to acurately 
predict deformation, there is still a question whether a macro-model can be used for failure and 
life prediction because these are known to depend on local events. However, with a micro-model 
we can examine local behavior. Consequently, we show the local equivalent inelastic strains 
corresponding to locations marked A-M on Fig. 1 lb and N-S on Fig. 1 la tabulated in Table 4. 
The local equivalent inelastic strains are observed to vary greatly with loading and location. 
Transverse compression (points J, K, and L, Fig. lib) creates the largest localized inelastic 
strains shown in Table 4; 47jx, 1280)1, and 2490(1 for 1, 20, and 35 Pa/sec, respectively. The 
largest local-to-macro inelastic strain ratio is 4.7 and occurs for the 1 Pa/sec SCDR under 
transverse compression. The minimum subcell inelastic strain for each surface was zero. For 
transverse tension (points D, E, and F in Fig. 1 lb and Table 4) the macro inelastic strain is quite 
large due to debonding, while the local inelastic strains are either zero or relatively small. 

It is known that the overall inelastic strain rate, and hence the dissipation rate, has a 
discontinuity associated with the initiation of debonding (as well as when the matrix threshold 
stress is attained), this discontinuity is shown in Fig. 12 for points E and M on the 20 Pa/sec 
SCDR. Also shown in Fig. 12 is the increase in the overall (macro) J 2 with the applied strain £ 22 , 
clearly the values of Jj are much larger at point M than at point E due to the presence of an axial 
stress. The concavity in the 20 Pa/sec SCDR at point E in Fig. 1 lb is associated with the 
redistribution of fiber stresses to the matrix upon debonding, as such it is an artifact of the type of 
interface model employed and consequently should not be alarming. 

Inelastic surfaces can be determined experimentally at elevated temperature by applying axial- 
torsional loading to a thin-walled tubular specimen. It is desirable to determine an entire 
inelastic surface using a single specimen to avoid specimen-to-specimen scatter in results and 
make the program economically feasible. This requires absolutely minimizing the amount of 
inelastic deformation that occurs in each probe of the surface as any theory would require that 
these tests be performed at a fixed inelastic state of the material. Tests have been successfully 
conducted for target values of up to 100|i of permanent strain in each probe direction for 
monolithic materials containing no internal damage (e.g., Battiste and Ball, 1988). 
Consequently, in composite testing we must be aware of both the expected macro and micro 
inelastic strain fields developed during a given probe. As a result, we have tabulated in Table 4 
the predicted local equivalent inelastic strains that accompany the SCDRs in the On-On stress 
space for points N-S denoted in Fig. 11a. For the applied strain rate of lOOp/sec it appears 
feasible to experimentally determine a 1 Pa/sec SCDR on a single specimen as both the local and 
macro equivalent inelastic strains should be relatively small at points N and Q. However, it may 
not be possible to experimentally determine a complete 1 Pa/sec SCDR on a single specimen 
because of a change in material state associated with overshooting the target value that would be 
likely to occur due to the relatively large (predicted) discontinuity in the dissipation rate (e.g., 10- 
15 Pa/sec shown in Fig. 12) associated with the initiation of debonding. Further, there is a wide 
disparity in macroscopic equivalent inelastic strain for shear and axial compressive loading given 
a value of dissipation rate. The ratio of macro equivalent inelastic strain between shear and axial 
compression for SCDRs of 1, 20, and 35 Pa/sec is 4.5 (N/Q), 9.5 (O/R), and 18.6 (P/S), 
respectively (Table 4). This wide range of equivalent inelastic strains could preclude the use of a 
single specimen to experimentally determine SCDRs, especially for relatively large dissipation 
rate target values. It should also be reaffirmed that the inelastic strain rate, and thus SCDRs and 
SCISRs, is very sensitive to the loading rate due to matrix viscoplasticity. 
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The shapes of the SCDRs in the three stress spaces (Fig. 1 1) indicate that axial-torsional tests 
on longitudinally reinforced thin-walled tubular specimens are not the most discriminating tests 
for validating the mathematical construction of the dissipation potential. Axial-transverse load 
probing experiments conducted on an in-plane biaxial test machine, (while being more difficult 
to conduct due to stress and strain gradients in the plane of the plate, coupling between axial and 
transverse strains, and the large potential for buckling in compression), would provide more 
discriminating and valuable information for the construction and validation of the mathematical 
representation of the dissipation potential. Biaxial transverse loading, G22-G33, would also 
provide valuable information, although specimens would be costly and difficult to manufacture. 
Alternatively, circumferentially-reinforced thin walled tubes could be employed to obtain 
surfaces in the O12-G22 stress space, and even though this stress space was not considered in the 
present work these surfaces are likely to have complex geometries because both shear and 
transverse loading result in debonding. 

Figure 13 shows 1,15 and 25 p/sec 2 SCISRs at 500°C, again the target values were chosen to 
provide the largest possible range of SCISRs. SCDRs and SCISRs are quite similar for small 
values, 1 Pa/sec and 1 p/sec 2 , respectively, as can be seen by comparing Fig. 11 and 13. The 
reason for the similarity appears to be that these small values really correspond to the initiation of 
either debonding or matrix inelasticity, causing the difference between and £jy£jy to be 

small. For larger target values, the shapes of SCDRs (Fig. 11) and SCISRs (Fig. 13) become 
quite different in all three stress spaces, particularly the G22-O33 stress space. It should be pointed 
out that since Zy-ejy and finely are different quantities, one should not expect a one-to-one 

correspondence between SCDRs and SCISRs. However, the values for SCDRs and SCISRs 
plotted in Fig. 1 1 and 13 were chosen to provide the widest possible range for the given loading 
rate and temperature, facilitating comparisons of trends between SCDRs and SCISRs for small, 
medium and large amounts of inelasticity. It is apparent that SCISRs are more affected by 
debonding than are SCDRs. This could be related to the discontinuity in the inelastic strain rate 
associated with the initiation of debonding. For a SCDR the discontinuous €y is multiplied by 

the effective stress Zy, while for a SCISR the discontinuity is magnified by multiplying e\j by 
itself. 

Local and macro equivalent inelastic strains associated with the SCISRs shown in Fig. 13 are 
presented in Table 5. Typically, both the local and macro equivalent inelastic strains are lower for 
SCISRs than for SCDRs. In some cases the macro equivalent inelastic strain is considerably 
smaller for SCISRs than for SCDRs (transverse tension and shear). However, even for SCISRs 
the local equivalent inelastic strains can be quite high in some cases (axial tension/compression), 
having significant implications for experimental programs. 

CONCLUSIONS 

In this study we have generated macro flow/damage surfaces using three distinct definitions 
for these surfaces, i.e., “yield” (deviation from proportionality), surfaces of constant inelastic 
strain rate (SCISRs), and surfaces of constant dissipation rate (SCDRs). These surfaces have 
been calculated using the method of cells micromechanics model for a representative silicon 
carbide/titanium metal matrix composite (i.e., SCS-6/TIMETAL 21S) at temperatures ranging 
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from 23 to 500°C. The micromechanical analysis explicitly accounts for matrix viscoplasticity, 
fiber-matrix debonding, and residual stresses. The main conclusions of this study are listed here. 

• The various macro flow/damage surfaces possess differences in tension and compression 
with their centers being offset from the origins as a result of the presence of internal 
(microscale) residual stress fields. 

• The shape of these surfaces are distinctly affected by the presence of internal damage (i.e., 
fiber-matrix debonding) as well as the applied biaxial stress state. 

• The size of surfaces of constant deviation from proportionality (SCDFPs) and SCISRs is 
decreased by debonding for transverse tension and shear loadings, distorting the shape of the 
inelastic surface. SCDRs are also affected by debonding, but to a lesser extent. 

• The range of possible values for SCDRs determined by various applied loadings is limited by 
the disparity in stress magnitudes at the onset of inelasticity. This disparity in stress 
magnitude is attributed to the presence of fiber induced anisotropy and fiber-matrix 
debonding. 

• Clearly axial/torsion tests are not the most discriminating tests for differentiating between 
surface definitions or illustrating the required mathematical constructs of the dissipation 
potential. Here, in-plane biaxial tests appear to be the more important and informative tests 
to conduct. 

• The magnitude of the pertinent macro target value significantly influences the magnitude of 
the resulting internal state (i.e., the magnitude of local inelastic strain and internal stress) and 
consequently the reusability of a single test specimen for multiple surface probes. 

• Given a macroscale starting viewpoint, relatively sophisticated mathematical representations 
(in the form of the invariants employed, the inclusion of Jj,, and possibly /j, etc.) for the 
macro dissipation potential would be required to capture the complex geometric properties of 
the various flow/damage surfaces. In contrast, employing a micromechanics model, such as 
the method of cells, one can use a relatively simple constituent constitutive (e.g., 7 2 
viscoplastic theory for the matrix) relation to obtain the required macro multiaxial behavior. 

• Fiber-matrix debonding is independent of axial loading, except possibly in the case of axial 
compression, where debonding could occur due to transverse Poisson expansion. 

• A key parameter in the ability to determine the entire SCDR accurately is the applied strain 
rate. 

The results of this study are in agreement with those from previous macro- and 
micromechanics studies relating to yield and flow surfaces for metal matrix composites 
exhibiting perfect bonding. Here we have included fiber-matrix debonding and paid particular 
attention to different target values as well as different definitions of overall flow/damage. 
Additionally, the micromechanics model has permitted us to compare overall and local stress and 
strain fields as well as provide guidance for an experimental program aimed at validating 
multiaxial constitutive relations for metal matrix composites. 
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TABLE 1: Ti-21s MATERIAL PARAMETERS 



23 °C 

300°C 

500°C 

650°C 

704°C 

E (GPa) 

114.1 

107.9 

95.1 

80.7 

59.7 

CTE (n°C') 

7.717 

9.209 

10.70 

12.13 

14.09 

K (MPa) 

1029 

768.4 

254.2 

5.861 

0.756 

(X (MPa/sec) 

667.6 

137.8 

1.45x1 O' 3 

6. 19x1 O' 9 

l.BxlO 11 

B 0 (MPa) 

6.908x1 O' 5 

1.035x1 O' 4 

2.756X10- 4 

5. 870x1 O' 4 

6.346x1 O' 4 

K (1/sec) 

0 

0 

1.68X10' 7 

1.00x1 O' 6 

6.01X10' 5 

P 

0.001 

0 

0 

0 

0 

temperature-independent: v = 0.365, 

n = 3.3, B\ = 

0.05, p = 1.8, q = 

1.35 


TABLE 2: INTERFACE PROPERTIES 


T 0 (MPa) 

u 0 (|im) 

T 1 

a 


200 

9 

1 

8 

0.4 
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TABLE 3: MATRIX EQUIVALENT INELASTIC STRAINS 
ASSOCIATED WITH SCDFPs SHOWN IN FIGURE 4 


Location 
(see Fig. 4) 

subcell 12 

(xlO -6 ) 

subcell 21 

(xlO -6 ) 

subcell 22 
(xlO" 6 ) 

maximum 

local/macro 

DFP 

A 

60.8 

17.7 

0 

6.1 

10(l 

B 

5890 

5350 

5190 

2.9 

0.2% 

C 

53.7 

0 

0 

5.4 

10(i 

D 

4860 

4460 

0 

2.4 

0.2% 

E 

40.7 

38.4 

18.6 

4.1 

10(i 

F 

5700 

5500 

5380 

2.8 

0.2% 

G 

0 

58.0 

0 

5.8 

10(1 

H 

2150 

7200 

0 

3.6 

0.2% 

I 

54.5 

0 

0 

5.4 

10(1 

J 

4620 

4700 

0 

2.4 

0.2% 

K 

41.6 

41.6 

0 

4.2 

10(1 

L 

5670 

5670 

5040 

2.8 

0.2% 

M 

41.4 

42.8 

0 

4.3 

10(1 

N 


5640 

5210 

2.8 

0.2% 

O 

0 

40.8 

0 

4.1 

10 (i 

P 

461 

6480 

387 

3.2 

0.2% 
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TABLE 4: MATRIX EQUIVALENT INELASTIC STRAINS 
ASSOCIATED WITH SCDRs SHOWN IN FIGURE 1 1 


Location 
(see Fig. 1 1 ) 

subcell 12 
(xlO -6 ) 

subcel!21 

(xlO -6 ) 

subcell22 

(xlO -6 ) 

macro DFP 
(xlO -6 ) 

maximum 

local/macro 

SCDR 

(Pa/sec) 

A 

24.6 

24.6 

17.7 

7.4 

3.3 

1 

B 

775 

775 

685 

246 

3.2 

20 

C 

1690 

1690 

1540 

542 

3.1 

35 

D 

0 

0 

0 

18.5 

0 

1 

E 

1.3 

0 

0 

1000 

0.001 

20 

F 

739 

0 

249 

4210 

0.2 

35 

G 

16.9 

16.9 

14.3 

5.1 

3.3 

1 

H 

618 

618 

574 

196 

3.2 

20 

I 

1350 

1350 

1270 

431 

3.1 

35 

J 

0 

47.0 

0 

10.0 

4.7 

1 

K 

149 

1280 

11.1 

307 

4.2 

20 

L 

504 

2490 

40.9 

641 

3.9 

35 

N 

0 

1.0 

0 

22.7 

0.04 

1 

O 

59.9 

17.4 

53.7 

1870 

0.03 

20 

P 

4140 

41.0 

4020 

8000 

0.5 

35 

Q 

16.9 

16.9 

14.3 

5.1 

3.3 

1 

R 

618 

618 

574 

196 

3.2 

20 

S 

1350 

1350 

1270 

431 

3.1 

35 
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TABLE 5: MATRIX EQUIVALENT INELASTIC STRAINS 
ASSOCIATED WITH SCISRs SHOWN IN FIGURE 13 


Location 
(see Fig. 1 3) 

subcell 12 

(xlCT 6 ) 

subcell 21 

(xlO -6 ) 

subcell 22 

(xlO -6 ) 

macro DFP 
(xlO"*) 

maximum 

Iocal/macro 

SCISR 

((i/sec 2 ) 

A 

14.1 

14.1 

9.6 

4.2 

3.3 

1 

B 

606 

606 

530 

192 

3.2 

15 

C 

1600 

1600 

1450 

511 

3.1 

25 

D 

0 

0 

0 

18.5 

0 

1 

E 

0 

0 

0 

18.6 

0 

15 

F 

0 

0 

0 

18.6 

0 

25 

G 

14.0 

14.0 

11.7 

4.2 

3.3 

1 

H 

649 

649 

602 

205 

3.2 

15 

I 

1740 

1740 

1640 

555 

3.1 

25 

J 

0 

14.3 

0 

3.0 

4.8 

1 

K 

28.4 

533 

1.6 

120 

4.4 

15 

L 

114 

1108 

8.2 

262 

4.2 

25 

N 

0 

1.0 

0 

22.7 

0.04 

1 

0 

0 

1.0 

0 

22.7 

0.04 

15 

P 

0 

1.0 

0 

22.7 

0.04 

25 

Q 

14.0 

14.0 

11.7 

4.2 

3.3 

1 

R 

649 

649 

602 

205 

3.2 

15 

S 

1730 

1730 

1640 

555 

3.1 

25 
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cr 22 (MPa) o 12 (MPa) o 22 (MPa) a l2 (MPa) 


SiC/Ti-2 1 s 
perfect bond 
no residual stresses 
overall flow surfaces 

micro-predicted, 1 Pa/sec SCDR 

macro-predicted, F = constant 





-2000 0 2000 -2000 0 2000 
a, [(MPa) a 33 (MPa) 

Figure 5: Comparison of micro- and macro-predicted overall 
flow surfaces at 23°C, no residual stresses 


SiC/Ti-2 1 s 
perfect bond 
residual stresses 
overall flow surfaces 

micro-predicted, 1 Pa/sec SCDR 

macro-predicted, F = constant 
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Figure 6: Comparison of micro- and macro-predicted overal 
flow surfaces at 23°C, residual stresses included 
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Figure 7: Effect of loading (strain) rate on SCDFPs and SCDRs 
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